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ABSTRACT 

Exploring the origin of Lya nebulae ( "blobs" ) at high redshift requires measurements of their gas 
kinematics that are impossible with only the resonant, optically-thick Lya line. To define gas motions 
relative to the systemic velocity of the blob, the Lya line must be compared with an optically-thin 
line like Ha A6563, which is not much altered by radiative transfer effects and is more concentrated 
about the galaxies embedded in the nebula's core. We obtain optical and n ear-infrared (NIR) spectra 
of the two brightest Lya blobs (CDFS-LAB01 and CDFS-LAB02) from the lYang etaTI pOlph sample 
using the Magellan/MagE optical and VLT/SINFONI NIR spectrographs. Both the Lya and Ha lines 
confirm that these blobs lie at the survey redshift, z ~ 2.3. Within each blob, we detect several Ha 
sources, which roughly correspond to galaxies seen in HST rest-frame UV images. The Ha detections 
show that these galaxies have large internal velocity dispersions (a v = 130 — 190 km s -1 ) and that, 
in the one system (LAB01), where we can reliably extract profiles for two Ha sources, their velocity 
difference is Av ~ 440 km s -1 . The presence of multiple galaxies within the blobs, and those galaxies' 
large velocity dispersions and large relative motion, is consistent with our previous finding that Lya 
blobs inhabit massive dark matter halos that will evolve into those typical of r ich clusters today and 
that the embedded galaxies may eventually become brightest cluster galaxies (jYang et al.ll2010r ). To 
determine whether the gas near the embedded galaxies is predominantly infalling or outflowing, we 
compare the Lya and Ha line centers, finding that Lya is not offset (Ai>Ly Q = +0kms~ 1 ) in LAB01 
and redshifted by only +230 km s -1 in LAB02. These offsets are small compared to those of Lyman 
break galaxies, which average +450kms _1 and extend to about +700kms~ 1 . In LAB02, we detect 
C II A1334 and Si II A1526 absorption lines, whose blueward shifts of ~ 200 kms -1 are consistent with 
the small outflow implied by the redward shift of Lya. We test and rule out the simplest infall models 
and those outflow models with super /hyper- winds, which require large outflow velocities. Because of 
the unknown geometry of the gas distribution and the possibility of multiple sources of Lya emission 
embedded in the blobs, a larger sample and more sophisticated models are required to test more 
complex or a wider range of infall and outflow scenarios. 

Subject headings: galaxies: formation — galaxies: high-redshift — intergalactic medium 



1. INTRODUCTION 

Lya nebulae, or "blobs," are extended sources at 
z ~ 2-6 with typical Lya sizes of > 5" (> 50 kpc) 
and line luminosities of Lj, vn > 10 43 ergs" 1 
(e.g.. Keel et al.l Il999t iSteidel et ail 120001: iFrancis et al.l 
20011 iMatsuda et al I T2004 120101 : iDev et al.l 12005 



Saito et al.l 120061: ISmith fc Jarvisl 120071: iHennawi et al 
2009 1 ; [Ouchi et *£j]|2009t IPrescott et al.ll2009l lYang et al 
20091 boiOf) . Their number density and its large variance 
suggests that blobs lie in massive (Mhaio ~ 10 13 M Q ) 
dark matter halos, which will evolve i nto those typical of 
rich clusters today (jYang et alj [2010). Galaxies embed- 
ded within blobs are likely to become brightest cluster 
galaxies. Therefore, blobs are important sites for study- 
ing the early interaction of galaxies with the surrounding 
intergalactic medium (IGM). 
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This interaction is probably tied on some scale 
to the source of the blobs' extended Lya emis- 
sion, but the mechanism is poorly understood. For 
example, the blob may represent emission from 
;alactic-scale outflows gen erated by star formation 
Taniguchi fe Shioval I2000D. inten se rad iative feed- 
back from AGNs (jHaiman fc Reesl 120011: IGeach et all 
2009), or even co oling radiation from accreting gas 
' Fardal et all 120011: IHaiman. Spaans. fc Quataertl I2000t 
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Diikstra fc Loebll2009l; iGoerdt et al.ll2010D . 

There are generally two ways to unravel the nature 
of the blobs. The first approach is to study what pow- 
ers the Lya emitting gas, counting all ionizing sources 
within or around the extended Lya emission and compar- 
ing that energy budget with the observed Lya luminosity. 
Sometimes, one can identify possible e nergy sources, e.g. , 
powerful AGNs d etected in X-rays (|Yang et al.l 120091 : 
IGeach et al.l 12009ft . Some potential sou rces may even 
lie outside a blob (|Hennawi et al.l I2009D . If no such 
source is found, the remaining Lya lu minosity may be at- 
tributable to cooling radiation (e.g.. iNilsson et alj|2006t 
ISmith fc Jarvisll200 7t ISmith et alj|2008l) . 

An alternative path is to observe spectroscopically the 
kinematics of the extended Lya-cmitting gas. Even ad- 
dressing whether the gas is falling into an embedded 
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galaxy or outflowing into the IGM would be a significant 
step forward, providing evidence for either gas accretion 
by forming galaxies or for AGN/starburst-driven galactic 
winds, respectively. Combining the two approaches by 
considering the blob's energetics and kinematics simulta- 
neously is ideal, a test of whether the bulk motions of the 
gas are in fact coupled to the source of its illumination. 
For example, ga s infall is possible even i n the presence of 
an AGN (e.g.. iWeidinger et all 120051: Humphrey et al.l 
120071: lAdams et al.l l2009f ). Cases in which blob gas is 
both outflowing on small scales due to galactic processes 
and infalling on larger scales from the extended IGM are 
not unlikely. 

Unfortunately, the kinematic approach has been 
stymied by the limitations of the Lya line, which, al- 
though bright compared with other lines, is hard to in- 
terpret. The ambiguity arises because Lya is a reso- 
nant line and typically optically thick in the surround- 
ing intergalactic medium. As a result, studies even of 
the same blob ' s kine matics can disagree. On one hand, 
IWilman et aD (2005) argue that IFU spectra from a cer- 
tain blob are consistent with a simple model where the 
Lya emission is absorbed by a foreground slab of neutral 
gas swept out by a ga lactic scale outflow. On the other, 
Diiks tra et al.l (|2006bO claim that the same data can be 
explained by the infall of the surroundi ng intergalactic 
medium. Worse, iVerhamme et al.l ()2006l ) comment that 
the same Lya profiles are most consistent with static 
surrounding gas. 

To distinguish among such possibilities requires a com- 
parison of the Lya line center with the c enter of an 
optic ally-thin line like HaA6563 (see also lYang et all 
2006). Ha is a better measure of the blob's systemic 
velocity, i.e., of the precise redshift, because it is not seri- 
ously altered by radiative transfer effects and is more con- 
centrated about the galaxies in the blob's core. Despite 
the complexity of various blob models, infall models pre- 
dict that the peak of the optically thick Lya emission line 
must be blueshifted with respect to Ha, wherea s outflows 
will redshift the Lya line (Diikst ra et al.ll2006al) . This di- 
agn ostic has been applie d both to local starburst galax- 
ies (|Legrand et al.l 119971) and to Lyman break galaxies 
(LBG s) at z =1.5-2.5 (jPettini et al.l [20011: iSteideleFall 
2004) and has provided evidence for outflows driven by 
strong stellar winds. An analysis in a similar spirit, but 
using the H I 21cm absorption line to mark the systemic 
redshift of a high-z radio galaxy, finds that the Lya line 
is blueshift ed and thus that ex tended gas is inflowing to 
the AGN (| Adams et al.1 120091 ) . Previous spectroscopic 
obser vations of Lya blobs (e.g.. at z 3: IWilman et al.l 
120051: iMatsuda et al.ll200llSaito et al.ll2008D include only 
Lya, precluding a direct comparison with a cleaner line's 
center and muddling the interpr etation of the kinematics . 

To overcome this challenge, lYang et al.l (|2009L l2010h 
conduct blind surveys for Lya blobs at z ~ 2.3, a redshift 
at which important rest-frame optical diagnostic lines 
(e.g., [O II]A3727, [O III] A5007, H/3A4868, HaA6563) 
fall in NIR windows and avoid bright OH sky lines (or at- 
mospheric absorption). In this paper, we present optical 
and NIR spectroscopy of the two brightest bl obs, CDFS- 
LAB01 and CDFS-LAB02, of the 25 in the lYang et al.l 
(2010) sample. For each, we compare the Lya line profile 
with the center o f the Ha line and use existing radiative 
transfer models (Diik stra et al.l l2006al IVerhamme et al.l 



120061 I2008T) to test whether the detected gas is infalling 
or outflowing. 

In addition to the blob's systemic velocity, NIR spec- 
troscopy can constrain (1) the velocity dispersion and 
any ordered motion like disk rotation of the brightest 
embedded galaxies, (2) the presence of AGNs with rest- 
frame optical nebular line diagnostics (e.g., [N II] /Ha), 
(3) the gas-phase metallicity, and (4) the star formation 
rate (SFR) directly from the Ha flux, allowing a compar- 
ison with that derived from the rest-frame UV. Before 
now, none of these measurements has been available for 
Lya blobs. In this paper, we focus on the kinematics of 
the gas in the blob and in the embedded galaxies, defer- 
ring discussions of the AGN line diagnostics, gas-phase 
metallicity, and SFRs to a future paper (Y. Yang et al. 
2011, in preparation). 

In §^J we review how we selected blobs fro m our nar- 
rowband imaging surveys (|Yang et all [2010) and then 
describe the optical/NIR spectroscopy of the first two 
blobs. In Sj3j we present the first results from the spec- 
troscopic campaign. We confirm spectroscopically the 
blob redshift and discover multiple Ha sources within 
each blob (33TTJ) . Wc compare the Lya profiles with the 
Ha line centers to discriminate between simple infall and 
outflow scenarios (i j3.2p . In ij3.3l we describe the inter- 
stellar absorption lines detected in LAB02 and compare 
the low ionization lines with Ha. We estimate the ve- 
locity dispersions of three of the embedded galaxies in 
§3.41 In §31 we compare the Lya profiles and Lya— Ha 
offsets with existing radiative transfer models. We also 
discuss our results in the conte xt of the clumpy circ um- 
galactic gas model proposed by iSteidel et al.l (|2010D . §3] 
summarizes our conclusions. Throughout this paper, we 
adopt cosmological parameters: Hq = 70kms _1 Mpc -1 , 
Q M = 0.3, and Q A = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Sample 

We observe the two brightest Lya blobs from the 
I Yang et"aTI (|2010T ) sample, which consists of 25 blobs 
discovered in four different survey fields. These two 
blobs lie in the Extended Chandra Deep Field South 
(E-CDFS) and were found via deep imaging with the 
CTIO-4m MOSAIC-II camera and a custom narrowband 
filter (JVB403). This filter has a central wavelength of 
A c w 4030A, designed for selecting Lya-emitting sources 
at z 2.3, and a band-width of AAfwhm ~ 45A, provid- 
ing a line of sight (hereafter LOS) depth of Az ~ 0.037 
that corresponds to 46.8 comoving Mpc. The blob selec- 
tion criteria are: EW b s > 100A and isophotal area A[ so 
> 10 □" above the surface brightness threshold of 5.5 x 
10~ 18 ergs -1 cm~ 2 arcsec -2 . We show the images of the 
two blobs (CDFS-LAB 01 and CDF S-LAB02) in Figure 
[TJ We refer readers to lYang et all (|2010l ) for details of 
the sample selection. 

2.2. Optical Spectroscopy 

We obtained high resolution optical spectra of the two 
Lya bl obs using the Magell an Echellette Spectrograph 
(MagE; lMarshall et al.l f2008) on the Magellan Clay 6.5m 
telescope. MagE provides a spectral resolution of R ~ 
4100 for a 1" slit over a wide wavelength range of 3000A 
- lOOOOA. While the slit length is small (10") compared 
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Fig. 1. — (Left) Images of Lya blobs at various wavelengths: U, continuum-subtracted Lya line, B, K, Spitzcr IRAC 3.6/^m, and HST 
F606W images. Contours represent a surface brightness of 4x 10 -18 ergs" 1 cm -2 arcsec - 2 from our CTIO-4m narrowband imaging. Ticks 
represent 10" (82 physical kpc) intervals. The red crosses in the IRAC 3.6/im image indicate X-ray sources from the CDFS and F-CDFS 
catalogs. (Right) High-resolution HST and SINFONI Ha images collapsed in the wavelength direction at A ~ 2.175/^m. Here, we show the 
SINFONI observation taken in the best seeing for LAB01 (1.3 hrs with 0"3 in K-band). For LAB02, the seeing was ~l'/2. SINFONI reveals 
two and three relatively bright Ha sources within LAB01 and LAB02, respectively. These sources roughly correspond to embedded galaxies 
or galaxy fragments seen in the HST images (labeled as A, B, C and D). The gray and green boxes represent the MagE slit positions and 
the extraction windows for the spectra shown in Figure [5] respectively. In the SINFONI image, we show the final FOV resulting from the 
dither patterns (dot-dashed lines). The negative images outside the FOV are the ghost images due to sky subtraction. 



to typical longslits, the excellent blue sensitivity around 
4000A is ideal for obtaining the Lya line profile. Obser- 
vations were carried out on UT 2008 July 26 and October 
1, and 2009 February 17. For LAB01, we first acquired a 
low resolution spectrum in 2008 July 26 to spectroscopi- 
cally confirm the redshift prior to the NIR spectroscopic 
run. Later, we obtained the R ~ 4100 spectrum with a 
slit orientation of P.A.=42° as shown in Figure [T] (Lya 
panel). For LAB02, we used a 2" slit and a P.A.=0° 
to obtain aJ?~ 2000 spectrum. The sky condition was 
clear, and the seeing ranged from 0.9 to 1.3". The indi- 
vidual exposure times varied from 20 to 40min, depend- 
ing on the variability of the seeing. Total exposure times 
were 2.5 and 2.0 hrs for LAB01 and LAB02, respectively. 
For accurate wavelength calibration, we took ThAr lamp 
frames right before and after the science exposures and 
at the same telescope pointing to eliminate systematic 
errors introduced by instrument flexure. 

We reduce the MagE data usi ng the Danie l Kel- 
son's Carne gie-Python Package (Kel son et al.l 120001 : 
IKelsonl [20031 ) . The frames are overscan-corrected, bias- 
subtracted, and flat-fielded with Xenon-flash and quartz 
(incandescent) lamps. Using night sky lines and ThAr 
lamp spectra, we estimate the uncertainties in the wave- 
length solution to be ~0.03A (0.1 pixel). The sky back- 
ground is then subtrac ted in 2D using the method de- 
scribed in lKelsonl (j2003), and the frames are corrected for 
spatial distortions. We combine the rectified 2D spectra 



and extract the ID spectra using IDL routines. We also 
re-reduce the data with the MagE S pectral Extractor 
package (MASE; IBochanski et al.ll2009fl and confirm that 
the results are consistent. The spectra are flux-calibrated 
with 2-3 standard stars each night. Due to the limited 
spatial coverage of the slit and the relatively short expo- 
sure times, our spectra are dominated by the brightest 
spots within the blobs, which coincide with the galaxies 
detected in the HST rest-frame UV images. Future IFU 
observations should be helpful in constraining how the 
Lya profiles vary over the blobs. 

2.3. Near Infrared Spectroscopy 

We carried out NIR Integral Field Unit (IFU) 
spectroscopy of the two Lya blobs usin g SINFONI 
(jEisenhauer et al.H2003t IBonnet et aT]|2004D on the VLT 
UT4 telescope in visitor mode between UT 2009 January 
30 and February 4. To measure the redshifted Ha (A b s 
~ 2.175/^m) line, we used the K-band grating, which cov- 
ers 1.95/im - 2.45/im with a spectral resolution of R ~ 
4000 (Av ~ 75kms _1 ) and dispersion of 2.5Apixel _1 . 
Because the blobs are spatially extended over 5"— 10" 
and include multiple sources, we employed the seeing- 
limited mode (i.e., no adaptive optics), which provided 
the largest FOV (8" x 8") and a plate scale of 0"25 pixel" 1 
for the spatial resolution elements (spaxels) . In each ob- 
serving block, we observed reference stars, which we used 
to blind offset to the blobs as well as to monitor the tele- 
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Fig. 2. — Integrated ID Lya spectra for LAB01 and LAB02 at z ~ 2.31. The spectra were extracted from apertures of l"x2"4 (LAB01) 
and 2"x2'.'4 (LAB02). (Right) MagE spectra in the observed wavelength frame from 3500A - 7000A. Vertical dot-dashed lines indicate 
the location of important UV emission/absorption lines. The gray bars indicate the masked region due to strong sky lines. The spectra 
are boxcar-smoothed with 5 spectral pixels. In addition to strong Lya emission at ~4030A, LAB01 shows C IV A1549 and He II A1640 
lines. In LAB02, UV continuu m from the embedde d galaxies, possibly galaxies A and B+C, is detected. For comparison, we show the 
composite LBG spectrum from Shaplcy ct al. (2003). (Left) Close-up of Lya profiles in the rest-frame. We convert the observed spectra 
into rest-frame wavelength using the redshifts determined from their He* lines (z = 2.31421 and 2.31928). 



scope pointing accuracy and the seeing variations. The 
K-band seeing ranged from 0'/3 to 1'.'2, and the sky con- 
dition was clear. To maximize the on-source integration 
time, we adopted an "on-source dithering" strategy in 
which the science targets were always kept within the 
FOV but at different detector positions. For LAB01, we 
used a three-point dithering scheme that moves among 
the pointings by ^4" to south or west, providing a usable 
FOV of two 4"x4" squares (see the SINFONI panel in 
Figure [T]). For LAB02, we employed a two-point dither 
pattern in the north-south direction, resulting in a 8" x4" 
effective FOV (dot-dashed line in Figure [1]). The total 
integration times for LAB01 and LAB02 are 6 hrs and 
0.5 hrs, respectively, which consist of 5 or 10 min indi- 
vidual exposures depending on the sky condition. For 
the telluric correction and flux calibration, we observed 
four B or G stars every night before or after the science 
targets with wide range of airmass. 

We reduce the SINFONI data using the ESO pipeline. 
Dark current and sky background are removed from each 
science frame by subtracting the pseudo-sky frame, which 
is constructed from the average of the two science frames 
bracketing each science frame. Then, we flat-field the 
data and correct for bad pixels. Because we will com- 
pare the velocity centers of the Lya and Ha lines, accu- 
rate wavelength calibration is critical. We compare the 
wavelength solution obtained from the OH sky lines in 
the science frames to that from the daytime arc lamps. 
The arc lamp wavelength calibration error arising from 
instrument flexure can be as large as 1.5 A, depending on 
the airmass (which ranges up to ^2.0). Therefore, we 
adopt the OH sky line wavelength calibration, whose er- 
ror is typically ~ O.lA (0.04 pixel) for each science frame. 
Using this wavelength solution, we constructed 3D data 



cubes from the 2D spectra after correcting spatial distor- 
tion. Because it is not possible to identify /centroid the 
sources in individual exposures, we align them according 
to the dither offsets within each observing block (OB). 
The spatial offsets between OBs are calculated from the 
PSF calibration frames and are typically small (~ C/l). 
When combining the data cubes, we adopt an iterative a- 
clipping algorithm to reject possible outliers and produce 
the "sigma cubes" that represent the standard deviation 
of the adopted pixels. Because no stellar continuum is 
detected, we subtract the median sky values from each 
wavelength plane to avoid the strong variation of global 
sky background. Finally, we flux-calibrate the cubes us- 
ing the broad-band magnitudes of standard stars. Atmo- 
spheric absorption is minimal at Ha due to our choice of 
survey redshift, so we do not apply a telluric correction. 

3. RESULTS 

3.1. Redshift Confirmation with Lya and Ha Lines 

We spectroscopically confirm the redshifts (z ~ 2.31 
- 2.32) of the two Lya blobs using both the Lya and 
Ha emission lines. Figure [1] shows the two Lya blobs 
at various wavelengths: 17, continuum-subtracted Lya 
line, B, K, Spitzer IRAC 3.6/zm, and HST F606W im- 
ages overlain with the Lya contours (left six panels). The 
archival ground-based images ( U, B, and K) , the Spitzer 
IRAC images, and the HST images are obtained from 
the Multiwavelength Survey by Yale-Chile (MUSYC; 
ICawiser et aI1l2006aD . the Spitzer IRAC / MUSYC Pub- 
lic L egacy in E-CDFS survey (SIMPLE; iDamen et afl 
I2010T ) . and the Galaxy Evoluti on from Morphology and 
SEDs (GEMS; iRix et I all l200l survey, respectively. In 
the right panels, we show the enlarged HST images and 
SINFONI Ha intensity maps that are collapsed in the 
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Fig. 3.— Integrated ID spectra of the Ha and [N II] AA6549,6583 lines extracted from a 2"x2" aperture for LAB01 and LAB02. Note 
that to maximize the S/N, we adopt different aperture sizes here than those used for Lya in Figure [2] The dashed and dot-dashed lines 
represent Gaussian fits to [N II] and Ha, respectively. The gray histograms in the lower panels show the uncertainties. For LAB01, two sets 
of Ha+[N II] lines from galaxies A and B are indicated. We use all available data (6 hrs) taken in various seeing (0''3 - l'/2). These two 
galaxies have a velocity difference Av (A-B) = 440 km s^ 1 . In our subsequent analyses, we compare only the Lya and Ha lines extracted 
along the line of sight to galaxy A. For LAB02, we are not able to reliably extract individual profiles from galaxies B and C, because of 
the low S/N. Thus we fit only one Ha+[N II] set to the LAB02 spectrum, which is dominated by galaxy A. 



wavelength direction centered at A ~ 2.1746/im (AA ~ 
66A) and 2.1788^ (AA ~ 42A) for the two blobs, re- 
spectively. For the SINFONI image of LAB01, we show 
only data taken in the best seeing (~ 0'.'3) to be compa- 
rable with the HST resolution. 

Interestingly, in the HST rest-frame UV images, 
both Lya blobs include multiple galaxies and/or frag- 
men ts, some of which appear to be interacting (see 
also iColbert et all I2006D . A few are also detected in 
Ha, indicating th at they lie at the same redshift. In 
I Yang et"ail (|2010l ) , we found that the number density and 
the field-to-field variance of Lya blobs are consistent with 
their occupying massive dark matter halos (~ 10 13 iW©) 
that probably evolve into those of typical galaxy clus- 
ters today. Thus, while some clum ps / fragments could 
be part of a one larg e galaxy (e.g., lLabbe et al.l 120031 : 
lElmegreen et al.l I2009D . we speculate that others may 
merge and evolve into present-day brightest cluster 
galaxies. 

Figure [2] shows the ID Lya spectra extracted from 
l"x2"4 and 2"x2"4 MagE slits for LAB01 and LAB02, 
respectively. In addition to strong Lya emission, LAB01 
also has C IV A 1549 and He II A 1640 emission lines, 
implying the presence of hard ionizin g sources (e.g, 
iPrescott et al.ll200alScarlata et al.ll2009f ). We defer anal- 
ysis of these lines and their relative strengths to a future 
paper (Y. Yang et al. 2011, in preparation), which fo- 
cuses on the multi-wavelength properties of the full blob 
survey and the implications for the sources of blob emis- 
sion. LAB02 has several interstellar absorption lines that 
we use in §3.31 to help constrain the kinematics of gas 
around the embedded galaxies. 

In contrast to the optically thick and spatially ex- 
tended Lya emission, we can use the optically thin, more 
spatially concentrated Ha emission to make an indepen- 
dent determination of the blob redshift from the galaxies 
embedded in the centers of the blobs. We extract the Ha 
spectra from a 2"x2" aperture in the SINFONI cubes. 
Figure [3] shows the resulting spatially integrated Ha pro- 



files within the blobs. The LAB01 Ha spectrum includes 
all 6 hr exposures combined to obtain higher S/N and to- 
tal flux. Note that whenever we show individual Lya or 
Ha spectra (Figures [5] and [3]) , we adopt different aper- 
ture sizes that maximize the S /N. On the other hand, we 
use the same aperture when comparing the Lya and Ha 
lines directly (§32]). 

For LAB01, we focus on the brighter SW component of 
the blob. The excellent seeing of our SINFONI observa- 
tion (~ 0.3") allows us to spatially resolve two embedded 
galaxies (A and B) in the SW part of the blob (hereafter 
LAB01A and LAB01B; Fig. [TjEl. These two galaxies are 
also kinematically separated as shown in the SINFONI 
Ha channel map (Figure[4]), a series of Ha intensity maps 
constructed by stepping through the 3D SINFONI data 
cubes in increments of 70kms~ 1 . Galaxy B shows an 
apparent velocity shear of ~ 250 kms -1 , suggesting disk 
rotation, while galaxy A has no obvious velocity struc- 
ture. In their spatially integrated Ha spectra, these two 
galaxies have the nearly same redshift, z£ a = 2.31421 ± 
0.0001 and z£ a = 2.30928 ± 0.0001, a separation of Az = 
0.0049 or Av(A-B) = 440 kms" 1 . This relative velocity 
is a poor estimate of the velocity dispersion a of the halo 
in which these galaxies are embedded. Nevertheless, it is 
consistent with the a expected for a massive, ^10 13 M Q 
halo. 

LAB02 also contains several continuum sources in the 
HST UV images (labeled as A-D). At least one of the 
brightest sources (LAB02A) is detected in the SINFONI 
Ha image. With SINFONI, we marginally detect the two 
sources (B+C) blended in the HST image. It is also pos- 
sible that galaxy D is blended with the much brighter 
galaxy A in the IFU data. We are not able to reliably 
extract Ha profiles from galaxies (B+C) and D because 

1 While it is possible that these two components could be H II 
regions or star-forming clumps in a single galaxy, it is unlikely given 
that their projected separation is ~5kpc and the velocity offset is 
~440 kms^ 1 . 
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Fig. 4. — Ha channel map for LAB01 and LAB02, respectively. The grids represent 1" intervals, and each channel is centered at the 
velocity noted at the upper left corner of each panel. Each intensity map is smoothed with a Gaussian kernel of FWHM=0'.'25. In LAB01, 
there are two components that are spatially and kinematically distinct, which we associate with galaxy A and B in the HST images. In 
LAB02, the Ha emission is dominated by galaxy A. There are two weak features that might correspond to galaxies B+C and D; because 
of the low S/N, we are not able to extract spectra for them. The last two channels of LAB02 are strongly affected by an OH sky line. 



of the shallow exposure (30 min) and insufficient spatial 
resolution. The integrated spectrum of LAB02 is domi- 
nated by the brightest galaxy A, for which we measure 
the zho = 2.31928 ± 0.0001. We summarize the proper- 
ties of Ha lines in Table [TJ 

While it is possible that the Ha detections include 
some extended gas, the centers of the emission coin- 
cide with the galaxies resolved in the HST image and 
are likely to originate there. Therefore, in the following 
section, we compare the Lya and Ha profiles along the 
LOS toward these embedded galaxies. There are small 
astrometry offsets (~0.2") between HSTUV images and 
SINFONI Ha maps (right panels in Figure [J). While it 
is likely that this mismatch arises from the astrometric 
calibrations of SINFONI, which rely entirely on blind 
offsets, we extract the Lya and Ha spectra from the 
same aperture, accepting the current astrometric solu- 
tion. Changing the relative position of extraction aper- 
ture (green boxes) for Lya and Ha does not affect any 
of the conclusions in this paper. 

We show the selected extraction apertures (small green 
boxes) in Figure [TJ In LAB01, we extract the spec- 
tra along the LOS toward galaxy A, because the peak 
of Lya map coincides with galaxy A. While galaxies A 
and B are well-resolved in the Ha image, the spectrum 
extracted from the sight-line to LAB01A (green box in 
Fig.Q} is slightly contaminated with light from galaxy B. 
Therefore, we simultaneously fit the profile with a total 
of six components (three for each galaxy) and remove the 
galaxy B components when comparing the Ha and Lya 
line centers in the following section. 

3.2. Comparison Between Lya and Ha Lines 

The Lya profile emerging from the blob should be 
complicated, a mixture of the bulk motions of the gas 
and geometry. For the details of Lya radia tive trans- 
fer, we refer readers to the literature (e.g.. lAhn et all 
200ll [20021: iDiikstra et alll2006aHVerhamme et al.l[2006l; 



Kollmeier et all 120101: iFaucher-Giguere et all 12010ft . In 
this section, we compare the Lya and Ha lines to con- 



strain the gas kinematics within the blobs. More specif- 
ically, we measure the velocity shift of the Lya line rel- 
ative to the Ha line center to distinguish between any 
infall and outflow of the IGM relative to the blob's sys- 
temic velocity. Once we measure this first-order kine- 
matic diagnostic (i.e., any blueshift or redshift of the Lya 
line), we discuss the detailed Lya and Ha profiles in Sj4j 
comparing them with radiative transfer calculations to 
determine the velocity of the dominant bulk motion. 

To illustrate what we expect from this experiment, we 
first consider the simplest geometric case: a galaxy sur- 
rounded by a spherically symme tric gas cloud that is 
static, collapsing or expanding ()Diikstra et all [2006a; 
iVerhamme et all [20061. Here we assume that the Ha 
line arises from the embedded galaxy, as we do not detect 
extended Ha emission, and that the Lya line originates 
either in the same region of the galaxy before eventu- 
ally escaping the surrounding cloud or in the cloud itself. 
Thus this blob will be observed as a point source in Ha, 
but as an extended source in Lya. Note that the Lya 
and Ha spectroscopy is along the LOS toward the cen- 
tral galaxy. While Ha from the galaxy can escape the 
surrounding gas without suffering absorption, Lya pho- 
tons from the galaxy and from gas along the LOS will be 
resonantly scattered many times until their frequencies 
are right for escaping the cloud, a condition dependent 
on the bulk motions of gas along the LOS. 

Because the cloud is optically thick to the Lya line, 
Lya photons on the Doppler wings should escape the 
cloud through random scatterings in the frequency do- 
main. If the surrounding gas is static, the Lya line will 
have a double-peaked profile, while the optically thin Ha 
photons will escape without any radiative transfer effects. 
If the cloud is collapsing, Lya photons within the cloud 
on the red side of double-peaked profile will see higher 
optical depth due to the LOS infalling gas, and the red 
peak will be depressed. If the gas is outflowing, the blue 
side of the profile will be more diminished. 

Therefore, if there is gas infall, the Lya profile will 
be asymmetric and blueshifted against the Ha line due 
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Fig. 5. — Comparison between optically thick Lya line (blue, solid) and optically thin Ha line (red, dashed) for the two Lya blobs. For 
both LAB01 and LAB02, we are comparing the Lyo profiles with Ha lines extracted from along the LOS toward galaxy A. Note that both 
Lya and Ha are extracted from the same aperture (green boxes in Figure [TJ for appropriate comparison. The Lya profiles are transformed 
to the observed NIR wavelength (SINFONI). For reference, we show the velocity frame obtained from the Ha redshift on the top x-axis. In 
both blobs, the Lya lines are broader than the Ha lines, implying that Lya is resonantly scattered by an optically thick medium. The Lya 
peaks either agree with or are redshifted against their Ha line center (A-UL ya = Okrns - 1 in LAB01 and 230kms — 1 in LAB02), suggesting 
no inflow and a possible outflow, respectively. 



to radiative transfer in the optically thick medium, but 
the corresponding Ha line profile should be symmetric. 
For outflowing gas, the Lya will be redshifted against the 
symmetric Ha line. The amount of shift depends on vari- 
ous parameters such as infall/outflow velocity and optical 
depth. Note that one cannot determine whether the Lya 
line is redshifted or blueshifted against the background 
velocity field unless there is a optically thin reference line, 
i.e., Ha. 

The actual geometry of the system is likely to be more 
complicated O- Yet the blue- or redshift of the Lya 
photons relative to the systemic velocity will occur as 
they make their way out of the system, as long as they en- 
counter opticall y-thick gas that is infa lling or outflowing 
along the LOS (jDiikstra et al.l l2006al ) . In other words, 
we are studying the bulk motion of gas along the pencil 
beam toward the embedded galaxies when we compare 
the Lya and Ha profiles extracted from this sight-line. 
In principle, it is even possible to compare the velocity 
offsets between extended Lya and extended Ha lines in 
two dimensions (i.e., Ha and Lya screens) with much 
deeper Ha spectroscopy and optical IFU data. 

In summary, the basic assumptions for our simple test 
using the velocity offset between the Lya and Ha lines 
are: 

1. Lya should go through an optically thick column 
of gas (either infalling or outflowing) to escape. 

2. Ha represents the systemic velocity of the entire 
system, and gas is moving (infalling or outflowing) 
relative to this center. 

With these assumptions in mind, we now consider the 
Lya and Ha line profiles for our two Lya blobs. We 
then examine these assumptions further and discuss how 
to overcome their limitations in []4j 

We show the Lya and Ha profiles for the two Lya 
blobs in Figure [5j We shift each blob's Lya profile into 
its Ha frame using Agf R = A^ cal x fffff. We place 
both blobs in the same velocity frame using the redshifts 



obtained from Ha. Note that all the optical and NIR 
wavelengths are converted to vacuum wavelengths and 
corrected to the heliocentric frame. In both blobs, the 
Lya line is broader than the Ha line and has a more com- 
plicated structure (i.e., an asymmetric profile or multiple 
peaks), implying that the Lya lines do experience reso- 
nant scattering. 

LAB01 has a Lya line width of Awfwhm = 
520 ± 40kms -1 whose peak coincides with that of the Ha 
line. While the blue side of the Lya profile (v < km s _1 ) 
agrees well with that of the Ha line, it has an extended 
red wing up to at least v ~ +600 km s -1 . It is not clear 
whether this red wing arises from the several poorly- 
resolved red peaks. LAB02 shows a double peaked Lya 
profile with a stronger re d peak, which is oft en observed 
in high-z Lya galaxies (|Tapken et all 120071) . LAB02's 
Lya profile is broad, extending from ^— 800 km s -1 to 
~+1000 kms^ 1 . The stronger red peak itself has a line 
width of Adfwhm ~ 420 kms -1 . LAB02's Lya profile 
also has a sharp absorption feature (or lack of emission) 
at v ~ Okms -1 that coincides with its Ha line center. 

To quantitatively measure the velocity offset between 
the Lya and Ha lines, we must first determine each line's 
center. For the asymmetric Lya line, we measure the 
wavelength at t he peak flux, which is cons istent with the 
methodology of ISteidel etaLl (<2004 [2010). Because the 
Lya profile is noisy, we smooth the spectrum with a box- 
car filter of 3 spectral pixels (Av ~ 40kms~ 1 ) and mea- 
sure the wavelength of the brightest flux. For the sym- 
metric Ha profiles in both blobs, we fit the spectrum with 
a Gaussian profile including the neighboring [N II] lines. 
We fit three Gaussian components with the same velocity 
width, but different intensities, centered at 6563 A (Ha), 
6549 A, and 6583A ([N II]). As mentioned earlier, due 
to the slight contamination of LABOlA's spectrum with 
light from LAB01B, we simultaneously fit the observed 
Ha profile for both galaxies and show only the spectrum 
of galaxy A in Figure 

The velocity offset, Av-^ ya is defined as the differ- 
ence between the Ha line center and the Lya peak 
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TABLE 1 
Properties of Ha Lines 



Target z L(Ua) cr(Ra) Au Lya 

(xl0 43 ergs _1 ) (kms -1 ) (kms -1 ) 

CDFS-LAB01 A 2.31421 ± 0.0001 1.16 ± 0.06 184 ±9 +0 ± 20 

B 2.30928 ± 0.0001 0.78 ± 0.05 134 ± 8 

CDFS-LAB02 A 2.31928 ± 0.0001 2.42 ± 0.11 152 ± 5 +230 ± 30 



wavelength. We find a negligible offset between Lya 
and Ha (Av hya ~ +0±20kms _1 ) for LAB01A, and an 
offset, At>Ly Q ^ +230 ±30 kms" 1 , for LAB02A. These 
smaller than those of LBGs, which range 
from + 150kms ~ 1 up to +900kms ~ 1 with an average 
of +450 kms" 1 (ISteidel et alJiMoh . Among this LBG 
sample, only 12% of galaxies have Av-^ ya smaller than 
225 kms" 1 , thus the probability of finding Av-^ ya < 
230 kms" 1 for two Lya blobs by chance is rather small 
(~1.5%) if our Lya-Ha offsets are drawn from the same 
AwLy Q distribution^. 

For our two blobs, Ai>L yct is either consistent with zero 
(LAB01) or somewhat redshifted (LAB02). Therefore, 
according to the simplest picture, we find no evidence 
for a blueshift and thus for gas infall in either blob. The 
absence of a redshifted AuL yQ in LAB01 does not neces- 
sarily exclude an outflow, because LABOlA's Lya profile 
has an extended red wing and/or multiple red peaks, 
which can be outflow signatures. If we interpret the ob- 
served AuLya as a proxy for outflow velocity, then, at 
face value, any outflows in LAB01 and LAB02 are weaker 
than those typical of LBGs. In gl we further consider the 
implications of our A«Lya results using radiative transfer 
calculations. 

3.3. Interstellar Absorption Lines 

In addition to the Lya emission line, the optical spec- 
trum of LAB02 shows several absorption lines (Figure 
[2]). These absorption lines provide an opportunity to 
constrain the kinematics of cold gas on the near side of 
the galaxy, i.e., between the galaxy and observer, and 
thus to verify the detection of the modest AuL yQ outflow 
signature in the previous section. Note that the systemic 
velocity from Ha spectroscopy is again critical. 

Because the S/N of the continuum in the LAB02 spec- 
trum is very low, we first re-bin the spectrum by four 
spectral pixels, which worsens the resolution of the final 
velocity bins to ~90 kms" 1 . Figure [6] shows the ve- 
locity profiles of six absorption lines whose peaks were 
detected at the > 3cr level in the rebinned spectrum. 
We show the absorption profiles of three low-ionization 
interstellar lines (Si II A1260.42, C II A1334.53, and 
Si II A1526.72) and three high-ionization lines (Si IV 
AA1393.76, 1402.77 and C IV A1548.20). For compari- 
son, we also plot the comp osite LBG spectrum (R ~ 800) 
from IShaplev et all f|2003T ) and the individual spectrum 
of MS1512-cB58 (R ~ 1300) which is one of the appar- 
ently brightest LBG whose light is ampl ified due to the 
gravitational lensing ()Pettini et al. 2000). Although the 
S/N is low, the velocity ranges spanned by LAB02's ab- 

2 We note that it is not clear whether the distribution of Avi jyct 
in Steidcl et al. (2010) is representative of all LBGs given its sharp 
break around 300 kms -1 . 



sorption lines roughly agree with these two templates. 

We find that all interstellar absorption lines are 
blueshifted against the systemic velocity, indicating that 
the intervening gas is outflowing, which is consistent 
with the observed AvLya outflow signature discussed in 
the last section. To estimate the velocity of the out- 
flowing material, we measure how much the centroids of 
these velocity profiles are offset from the Ha line: Avis- 
While the high-ionization lines (Si IV AA1393, 1402 and 
C IV A1548) have the highest S/N, they are known to be 
contaminated with broader absorption features arising 
from stellar winds from massive stars (jLeitherer et al.l 
[19951 IShaplev et al.ll200l . Therefore, we measure the 
centroid velocities (Avis) using the three low-ionization 
lines. Two low-ionization lines (C II A1334 and Si II 

A1526) have Au IS lOOkms" 1 to -2Q0kms _1 . (It 

is difficult to reliably measure the centroid of the Si II 
A1260 line.) The absorption profile of C II A1334 might 
extend up to |i> ma x| ~ 600 kms -1 blueward, although its 
low S/N precludes a reliable measurement. If we adopt 
its centroid velocity as an outflow velocity, then it is sim- 
ilar to that of LBGs, which average (wis) = — 164kms _1 , 
but can be as large as -500 kms" 1 (|Steidel et al.ll2010D . 

3.4. Velocity Dispersions of Embedded Galaxies 

By fitting Gaussian profiles to the integrated Ha pro- 
files (Figure [3]), we obtain the velocity dispersions, ana 
= 184±9kms" 1 and 134±8kms" 1 for the LAB01A 
and LAB01B, respectively. For LAB02A, <j H a = 
152 ±5 kms" 1 . These velocity dispersions are corrected 
for the instrumental resolution (er rcs ~ 32 kms" 1 ). 

The internal velocity dispersions of the embedded 
galaxies (a v = 130 - 190 kms" 1 ) are high compared 
to tho se of star- forming galaxies at z ~ 2. lErb et all 
(2006) find that LBGs at z ~ 2 have an average of (er) 
= 108 kms" 1 with a standard deviation of 41 kms" 1 , 
excluding galaxies with AGNs. Our lowest ohc* ~ 
130 kms" 1 corresponds to the upper ~30% percentile of 
their d istribution . In pa rticular, there is only one galaxy 
in the lErb et al~l (|2006l ) sample that has higher a than 
LAB01A does (a Ra = 190 kms" 1 ). 

4. DISCUSSION 

In this section, we compare the observed Lya profiles 
and their offsets from the Ha lines with th e predictions of 
Lya radiative transfer (RT) calcula tions ()Diikstra et al.l 
I2006al iVerhamme eTaLI [20061 [20081 hereafter D06, V06 
and V08, respectively). To test the scenario in which gas 
in the surrounding cloud cools and flows toward the blob 
center, we compare our data with the spherical collapse 
model in D06. For the picture in which gas outflows from 
the embedded galaxies in the blob's core, we employ the 
expanding shell model investigated by V06. Note that 
these RT calculations assume that infall and outflow oc- 
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Fig. 6.— Velocity profiles of absorption lines in LAB02. (Top) Low-ionization lines : Si II A1260, C II A1334, and Si II A1526. (Bottom) 
High-ionization lines : Si IV AA1393, 1402 and C IV A1548. The gray and blue histograms a re the unbinned and re-binned (by 4 pixels) 
histogram s, respectively. Red and magenta lines represent the composite LBG spectrum (Shaplcy ct al. 2003) and the MS1512-cB58 
spectrum (Pcttini ct al. 2001]), respectively. The dashed horizontal lines indicate the 3xrms noise of the continuum. We use the low- 
ionization interstellar (IS) lines to measure velocity offsets, Adjs > 
outflow velocity inferred from the Lya and Ha line comparison. 

cur in spherical symmetry, in a collapsing sphere for the 
inflow model and in an expanding shell for the outflow 
model. As we argued previously, even such simple models 
should reproduce the general behavior of the Lya pro- 
file as long as the basic assumptions listed in §3.21 are 
valid. In §4.3[ we examine those assumptions in detail 
and discuss the future work required to overcome these 
limitations. 

4.1. Absence of Blueshift: No Gas Infall? 



In 33.21 we found that the Lya lines of the two blobs are 
not blueshifted relative to the Ha line centers within the 
measurement uncertainties. The absence of a blueshift 
rules out the simplest infall model where a spherical gas 
cloud is collapsing onto a central Lya source or where 
the collapsi ng cloud itself is e mittin g Ly a photons. In 
this mo del, iVerhamme et al.l (|2006f ) and iDiikstra et al.l 
(|2006al) study Lya radiative transfer, which depends 
on various parameters such as the location of the Lya 
source, the infall velocity field, the column density of 
neutral hydrogen (JVhi), and the Doppler parameter b 
representing the thermal and turbulent motion of the 
gas (see Fig. 4 and 7 in D06 or Fig. 5-7 in V06). The 
key prediction of these RT calculations is that the Lya 
spectrum is double peaked with an enhanced blue peak, 
producing an effective blueshift of the Lya profile. While 
the blue peak might be suppressed by the intervening 
IGM after the Lya ph otons escape the collapsing cloud 
(jDijkstra et al.l l2006b| ) . this effect is unlikely to be sig- 
nificant here because the mean transmission blueward of 
Lya should be relative ly high at z ~ 2.3: (e~ T ) ~ 0.8 
(jMcDonald et a l. 2001), where r represents Lya optical 
depth. 



The blueshift of Lya photons against the systemic ve- 
locity should arise as long as 1) there are infall motions 
along the LOS and 2) Lya photons on their way out of 
the system encounter optically thick gas with a velocity 
opposi te to their propagation direction (D iikstra et al.l 
2006a) . One possible way to hide a blueshift in this sim- 
ple infall model would be if the infalling gas is made 
optically thin by an ionizing source or has low enough 
column density so that the Lya photons do not suffer 
resonant scattering (see Fig. 11 in D06 or Fig. 6 in V06). 
In this case, we should still see nearly symmetric Lya 
profiles or extended blue wings, neither of which is ob- 
served in our Lya profiles. Therefore, we rule out this 
simple gas infall scenario for our two Lya blobs. 

In contrast to this simple model, the actual situation is 
likely to be more complex and the bulk motion of gas may 
not be spherically symmetric. For example, the Ha de- 
tections suggest that the galaxies embedded in the blobs 
are forming stars, which could generate mechanical feed- 
back into the surrounding gas cloud. There is also ample 
evidence for galacti c scale outflows in s tar-forming galax- 
ies at z = 2-3 fe.g-. ISteidel etalll2010h . Thus, one could 
imagine that even if gas infall takes place over larger 
radii (up to ^50 kpc), the innermost part of the gas 
cloud, close to the galaxies, could be strongly affected by 
outflows similar to those observed in other star-forming 
galaxies at z = 2 — 3. In this case, the emerging Lya 
profile will be more sensitive to the core of the blob, pre- 
sumably the densest part of the IGM, than to the gas 
infall. As a result, it might be difficult to detect the in- 
falling gas by measuring the Lya line shift. More simu- 
lations are required to test whether generic galactic scale 
outflows could affect gas infall and whether infall and 
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outflows can coexist. This kind of more realistic model 
has not been considered yet in RT calculations, so its 
spectral signatures (if any) are unknown. 

The other caveat is t hat the inflow might occur 
along filamentary streams (jKeres et al . 2005; Dcke l et al.1 

2009) and their covering fa ctor could be small 
(jFaucher-Giguere fc Kerei [201(1 such that we might 
completely miss the infall signatures in a particular LOS. 
We will further discuss this caveat in 34.3.11 

4.2. Small AvL ya : Weak or No Outflow? 

In §3.2| we found that the offset of the Lya line rel- 
ative to Ha was either consistent with zero (LAB01) or 
redshifted (LAB02). The existence of a modest outflow 
in LAB02 was further supported by the low-ionization 
absorption line shifts discussed in Section 1331 The ab- 
sence of a redshifted AuL yQ in LAB01 does not necessar- 
ily exclude an outflow because of its extended Lya profile 
toward the red. Here we compare the L ya profiles of our 
blobs with the RT outflow models of IVcrhamm e et al.l 
(|2006l |2008|) to interpret the velocity offset A^Lya in a 
more sophisticated manner. 

In these radiative transfer calculations, the Lya-Ha 
velocity offset is modulated entirely by the radiative 
transfer of Lya photons through a continuous medium 
within the shell. Another possibility, which we discuss 
in Section 14.3.31 is that the gas around the galaxy is 
clumpy and that the bulk motions of thes e clumps pri- 
maril y modulate the emerging Lya profile (jSteidel et al.l 

2010) . In the following section, we constrain the outflow 
velocity or "wind velocity" from the embedded galax- 
ies using the RT calculations. It is then possible to ask 
whether such an outflow is energetic enough (i.e., the ex- 
pansion velocity t> cxp is large) to expel/launch gas out 
to larger radii, thus creating the extended Lya emission 
observed as a blob with a typical size of > 50kpc. 

4.2.1. Comparison with Radiative Transfer Models 

To connect AvL yQ with an outflow velocity, we com- 
pare each blob's Lya profile with a model with an ex- 
panding shell geometry. Such a geometry would arise 
in the simple "superwi nd" or "hyperwind" mode l for 
Lya blobs prop o sed b y ITaniguchi fe Shioval (|2000H and 
iTaniguchi etatl (|2001| ). In this scenario, supernova ex- 
plosions or stellar winds following an intense starburst 
in galaxies develop into a so-called superbubble. If the 
kinetic energy deposited into the surrounding gas over- 
comes the gravitational potential energy of the galaxy, 
the gas clouds are blown out into intergalactic space as 
a superwind (e.g., Heckman et al. 1990). The gas blown 
out from the galaxies forms a shell, and the Lya photons 
are scattered from this expanding shell (with i> cxp ), es- 
caping preferentially by the scattering from the backside. 
In contrast, the Ha lines provide the systemic velocity of 
the star forming region. 

V06 and V08 demonstrate that the radiative transfer of 
Lya in this simple shell geometry is able to explain a wide 
range of observed Lya profiles. V08 classified Lya pro- 
files into three categories: (1) single-peaked asymmetric 
profile with an extended red wing, (2) profile with dou- 
ble peaks separated by sharp absorption at v ~ Okms -1 , 
and (3) asymmetric profile with bluer bump(s). Thanks 
to the Ha line, which provides the velocity center, we can 



determine that LAB01 and LAB02 correspond to (1) and 
(2) without ambiguity. 

We start with LAB01, which has an asymmetric Lya 
profile without any obvious blue peaks, but with an ex- 
tended red wing and/or multiple red peaks (i.e., category 
1 above). V06 show that a blue peak (v < Okms -1 ) is 
strongly suppressed with increasing (v exp /b) ratio, i.e., 
when the bulk outflow motion becomes more important 
than the thermal and turbulent motion of the shell gas. 
The resonant scattering in the expanding shell produces 
mainly two red peaks depending on the paths of the 
emergent photons: one peak near the line center (v ~ 
Okms -1 ) and another at v ~ 2u oxp , corresponding to 
photons that experience zero or one back-scattering from 
the receding shell, respectively. Here, the back-scattering 
of a photon means that it travels across the empty in- 
terior before re-entering the shell at a different location. 
Thus the location of the second red peak is a useful di- 
agnostic for the outflow velocity. 

Because of the degeneracy between many model pa- 
rameters (w eX p, b, iVui, idust), constraining the outflow 
velocity well requires fitting the spectra directly to the 
RT calculations. We do not have the means to do that 
here, so we estimate u oxp for LAB01 from how much its 
red Lya peaks are redshifted against its Ha line center. 
If A/hi is sufficiently large (Nm> 10 20 cm -1 ), the second 
red peak becomes larger and merges with the first red 
peak, so that the stronger red peak measures 2 v exp (see 
Fig. 1 of V08). If this case applies here, our AvL ya ~ 
Okms -1 implies w C xp^ Okms -1 , i.e., that the shell is 
static. 

On the other hand, decreasing -ZVhi, increasing dust 
extinction (rd us t), and decreasing v cxp /b tends to 
strengthen the first red peak, so that the brighter peak 
is located at < v < v cxp and does not carry any u cxp 
information. In this case, the location of the weaker sec- 
ond peak provides an estimate of t> exp . As discussed in 
£13.11 the LAB01 profile's extended red wing could be 
interpreted as multiple components (see Figure [SJ. If 
we adopt the small bump at v ~ 300 kins -1 as the sec- 
ond peak, the expansion velocity of the shell is v cxp ~ 
150 kins -1 . Considering both cases, the outflow velocity 
in LAB 01 could range from Okms -1 to ~ 150 kins -1 . 

The Lya profile of LAB02 is most consistent with 
the double-peak e d pro file (i.e, category 2 above) of 
iVerhamme et al.l |2008), who claim that such profiles 
arise when the shell is almost static (v cxp ~ Okms -1 ). 
Due to the Ha measurement, there is no ambiguity about 
the location of the velocity center. Even if we consider 
only the red side of the profile and classify it as a single 
asymmetric profile like LAB01, we would obtain an out- 
flow velocity only as high as ~ AwL yQ /2 = 115 km s -1 . 
The blueshifted centroids of the low-ionization absorp- 
tion lines in LAB02 are also consistent with this inferred 
outflow velocity. 

Therefore, if the expanding shell geometry is appli- 
cable to our two Lya blobs, we conclude that any 
outflows from the embedded galaxies are not strong 
{v C xp < 150kms -1 ) compared to the values (~ 500- 
1000 km s -1 ; see Eq. (2) of ITaniguchi fc Shioval l2000h 
required b y the extreme "supe r wind" or "hyperwind" 
scenarios (ITaniguchi fc Shioval 120001: ITaniguchi et~aT1 
120011 : lOhvama et al.1 120031 ) to explain the large size 
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(r > 50kpc) and broad Lya linewidth (FWHM ~ 
1000 km s" 1 ) of blobs. 

The simple shell geometry assumed here has some 
physical motivation and has been explored by exten- 
sive RT calculations (e.g., V06 and V08). However, this 
model implies that the outflows are intermittent with a 
small duty cycle, contrary to the observation that out- 
flows are seen in most star forming galaxies at high red- 
shift. We note that if there is a continuous stream of 
outflowing material, it might lead to a different interpre- 
tation of kinematic signatures discussed above. 

One other claim of an outflow in a blob was reported 
bv IWilman etaD pOOSh using IFU observations of the 
Lya line in Steidel's Blob 2. Using an expanding shell 
model, they found that the central wavelength, i.e., ve- 
locity, of absorption systems on top of the underlying 
broad Lya profiles is coherent across the whole blob with 
a blueshift of ^250 kms -1 . Because the systemic ve- 
locity of the blob gas is unknown in their case, there 
are other possible i nterpretations of their d ata, even in- 
cludin g an inflow (|Diikstra et al.l l2006bf) . IBower et al.l 
(|2004h and iWeiimans et al.l ([20091 ) found that Steidel's 
Blob 1 has an extremely complicated velocity structure, 
including velocity shears around the embedded galax- 
ies that can be interpreted as either outflow or rota- 
tion. It is interesting that all investigations of blob 
kinematics to date (in our CDFS-LAB01 and CDFS- 
LAB02, Steidel's Blobs 1 and 2) suggest relatively small 
outflow velocities. Therefore, the large wind speeds re- 
quired by super/hyper-wind models to explain the ex- 
tended Lya emission are excluded. One way to reconcile 
the small outflows with the "superwind" model is to as- 
sume that the outflow was once more energetic and that 
the expanding shell has been brok en up or slowed by 
contact with the surr ounding IGM (jWilman et al J 1200a 
IWeiimans et al.ll2009l ). Clearly, constraining the kine- 
matics of more Lya blobs is needed to address whether 
these small outflows are representative. 

4.2.2. Summary 

We conclude that the gas motion around the galax- 
ies embedded in these Lya blobs is not consis- 
tent with models like th e simple superwind picture 
(jTaniguchi fc Shiovdl200"oT ) that assume a spherical shell 
geometry and a shell velocity exceeding ^500 kms . 
Yet galactic scale outflows from star-forming galaxies at 
z = 2 — 3 appear to be common. The outflow properties 
have been extensively studied by comparing the Lya ve- 
locity center with those of optically-thin nebular lines 
and/or interstellar absorption lines (e.g.. iPettini et all 
l200lHShaplev et alJl2003HSteidel et al.ll2004L 120101 ) . Re- 
cently, iSteidel et al.l (|2010T ) show that Lya lines from 
LBGs are redshifted against Ha lines up to several hun- 
dreds kms -1 with an average of (AvL ya ) = +445 km s -1 , 
while the interstellar absorption lines are blueshiftcd 
against the nebular lines by (Auig) = —164 kms -1 . 

While these Av-^ ya and Avis offsets are interpreted 
as an outflow signature, it is not clear how the ab- 
sorbing material is distributed around galaxies. Were 
the gas blown out in a spherical sh ell, as in the super- 
wind picture discussed above (e.g. , iTaniguchi fc Shioval 
120051 ISchaerer fc Verhanmii 120081 V08), we could di- 
rectly compare the observed AvL ya of blobs and LBGs 



as a proxy for t> cxp . The velocity offsets of our two Lya 
blobs, A«Lya ~ Okms" 1 and 230 km s -1 , are small in re- 
lation to LBGs (|Steidel et al.ll2010l ): of 41 LBGs at z = 
2 - 2.6 with both Lya and Ha spectroscopy, ~88% have 
AwLya > 225 kms -1 . Therefore, any wind arising from 
galactic star formation has a lower expansion velocity in 
these two LABs than is typical of LBGs. 



4.3. Caveats 

We have adopted two basic assumptions in our analy- 
sis: that there is inflowing or outflowing gas in the LOS 
to the blob and that the Ha line marks the kinematic 
center of the blob gas. Here we examine how these as- 
sumptions may affect our co nclusions and comm ent on 
the recent model suggested bv ISteidel et al.1 (|2010f ) to ex- 
plain the circum-galactic kinematics of LBGs. 

4.3.1. Misalignment of Gas Flow to Line of Sight 

We have assumed that there is gas inflowing or out- 
flowing along the LOS to the embedded galaxies that 
lie in our spectroscopic slit. In this case, the observed 
Lya photons resonantly scatter in this optically thick in- 
falling/outflowing material until escaping into the LOS of 
the observer. However, we do not know how the blob gas 
is distributed around the embedded galaxies nor whether 
it is clumpy or smooth. 

For example, the ga s flow may not be iso tropic. As in 
bipolar outflows (e.g.. [Bland fc Tuflvl I1988D . a galactic- 
scale outflow may occur in the direction of minimum ISM 
pressure, which is typically perpendicular to the stellar 
disk. Or, the covering factor of the outflow material, the 
fraction of emission intercepted by the absorbers, could 
be rather high, at least 0.5 or up to u nity as observed 
in nearby AGNs dCrenshaw et al.H2003T) . I f gas accretion 
is taking place instead ( Keres et al.ll2005l l2009f ) , numer- 
ical simulations suggest that the infa ll may occur pref- 
erentially along fil amentary streams (jDekel et al.l 120091 : 
IGoerdt et afll2010l) . 

If the bulk motion of gas (cither infalling or outflow- 
ing) happens to be mis-aligned with our LOS, we may 
underestimate or even fail to detect the relative velocity 
shifts. Therefore, it is critical to expand the Lya-Ha 
comparison to a larger sample to average over any geo- 
metric effects and obtain better constraints on the inci- 
dence, direction, speed, and isotropy of bulk gas motions 
in blobs. If no large values of v exp are measured, even for 
a larger sample, outflows are unlikely to push the blob 
gas to larger radii, where it can be illuminated. On the 
other hand, an absence of blueshifted Lya lines would 
argue that infall is not dominant, at least not near the 
galaxies. 

Recent cos mological simulations includ ing radiative 
transfer (e.g.. iFaucher-Giguere et afll2010l ) are now able 
to examine the statistics of blobs whose extended Lya 
emission arises from cooling radiation. These models 
can predict the range of gas accretion covering factors, 
surface brightness profiles and ratio of blue-to-redshifted 
Lya lines for blobs. A larger sample of blobs with mea- 
sured Lya and Ha kinematics would provide a new test 
of these blob emission models. 
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4.3.2. Mistdentification of Kinematic Center 

Our second key assumption is that the Ha line marks 
the kinematic center of the surrounding gas, i.e., the ter- 
minus of any infall or the origin of any outflow. For an 
outflow, this assumption is reasonable because all pos- 
sible causes of the outflow (e.g., stellar winds, SNe, or 
AGNs) lie in or about galaxies, as does the Ha line it- 
self. For similar reasons, this assumption is likely to be 
valid for gas accretion (infall) into an embedded galaxy. 
However, it is possible that the accreting gas responds 
to the overall DM potential rather than to individual 
galaxies. Thus one could expect that detecting gas infall 
is more difficult than finding outflows using the technique 
presented in this paper. 

To overcome this issue, it is critical to increase the sam- 
ple size and measure how each blob's Lya profile varies 
spatially. This experiment should be conducted along 
multiple sightlines using using different Ha sources found 
in the same blob. In principle, it is possible to compare 
the offsets between the extended Lya and and the ex- 
tended Ha lines in two dimensions, i.e., between the Ha 
and Lya screens. Unfortunately, the Ha surface bright- 
ness is expected to be ~3% - 10% of Lya depending on 
the energy source. Thus, it will be challenging to obtain 
2D Ha IFU measurements in the near future. 

Currently, our Lya-Ha comparison includes only a 
handful of Lya blobs because of the limited redshift 
range where NIR spectroscopy of the optically-thin neb- 
ular lines (e.g., Ha and the [O III] A5007 line) is allowed. 
However, future CO observations with ALMA will also 
provide the systemic velocity of the embedded galaxies, 
permitting this test to be expanded to a larger sample of 
Lya blobs. 

4.3.3. Relationship to Clumpy, Circum- Galactic Medium of 

LBGs? 

As we were writing this paper, ISteidel et al.1 IpiM 
1201 it ) proposed a model for LBGs in which cool gas is 
distributed around the galaxy and accelerates radially 
outward with increasing galactocentric distance. They 
suggest that this blown-out material extends to ~80kpc, 
comparable to the typical size of Lya blobs, and that 
all LBGs would be classified as blobs if narrowband im- 
ages were sensitive to surface brightness thresholds 10 
times lower than typical of blob/emitter surveys (i.e., 
^10~ 19 ergs -1 cm -2 arcsec -2 ). Is this "clumpy circum- 
galactic medium (CGM) model" consistent with the 
properties of our Lya blobs? Unlike in the RT models 
that we considered earlier, the gas in this model is com- 
posed of small clumps (e.g.. lNeufeldlll99iT ; iHansen "Fog 
[200l . 

A critical part of the ISteidel et al.l (|2010f ) analysis in- 
volves stacking spectra of apparent galaxy-galaxy pairs 
at z ~ 2 — 3 to obtain a high S/N interstellar absorption 
profile arising from the gas surrounding the foreground 
galaxies. Sightlines to the background galaxies provide 
information on the spatial distribution of circum-galactic 
gas as a function of the galactocentric distances. There 
are no background galaxies behind our two blobs, and 
even constraining the absorption profile in LAB02A was 
challenging. We do not detect a stellar continuum in 
LAB01A or LAB01B. As a resul t, we are unable to test 
whether the ISteidel et ail (|2010|) LBG model could ex- 



plain the kinematics of either of our Lya blobs. Never- 
theless, we add a few q ualita tive comments here. 

In the Stcid el et al.l (|2010l ) model, Lya photons from 
the central galaxy are scattered off the surface of these 
discrete clumps; thus photons acquire a Doppler shift 
corresponding to the velocity of each clump. Lya pho- 
tons can escape the system when they achieve velocity 
offsets large enough to take the photons off-resonance for 
any material between the last scattering clump and the 
observer. Therefore, the Lya line shift is due to bulk 
motions of the clumpy material rather than to radiative 
transfer effects where Lya photons experience numerous 
resonant scatterings before escaping the cloud (Section 
14.21) . The profiles of Lya and the low-ionization inter- 
stellar absorption lines are determined by the velocity 
field and covering fraction of these gas clumps. The 
apparent peak of Lya emission is modulated primarily 
by the clumps with v ~ 0: Lya emission is more red- 
shifted and weaker for a wider velocity range spanned by 
the interstellar absorption. The velocity of the clumpy 
medium continually increases outward, unlike in the ra- 
diative transfer models discussed ea rlier, where the out- 
flowing shell had one velocity. In the ISteidel et al.l ([20100 
model, the velocity range spanned by the absorption lines 
reflects that of the surrounding gas distribution, and the 
Lya-Ha offset does not represent the outflow velocity. 

If this clumpy CGM model is applicable to blobs, then 
our observations of a small Lya-Ha offset could indicate 
an absence of neutral gas at v ~ rather than the weak 
gas outflow suggested by the RT models. One test would 
be to constrain the interstellar absorption line profile, 
which should extend to larger blueward velocities, but 
lack a component around v ~ Okms , if the clumpy 
model is accurate. Clearly, higher S/N spectra to reliably 
measure these weak absorption lines are required. 

The most notable outcome of the large spatial extent 
of a clumpy CGM around star-forming galaxies is dif- 
fuse Lya-emitting halos. Extended Lya emission has 
be en detected in stacked narrowband images of LBGs 
by IHavashino et al.l (|2004l) . More recently , by stacking 
deepe r narrowband Lya images of LBGs, ISteidel et al.l 
(| 2 1 If) show that all classes of star-forming galaxies ex- 
hibit diffuse Lya halos out to projected radii of ^80kpc, 
regardless of their spectroscopic properties (i.e., either 
absorption or emission in Lya within a spectroscopic 
slit), although these halos are 10 x lower surface bright- 
ness than typical blobs. The Lya surface brightness pro- 
files of the stacked images are similar among LBG sub- 
classes, with only different intensity scalings. Thus, Lya 
blobs may represent the most extreme example of galax- 
ies with diffuse Lya halos. The key difference of the 
clumpy CGM model is that the diffuse Lya halo origi- 
nates from Lya photons scattered from the galaxy's H II 
region(s) rather than from external mechanisms such as 
photoionization by AGN, photoionzation by metagalac- 
tic UV background radiation, or via cooling radiation. 

At this time, the applicability of the clumpy model to 
Lya blobs is unclear. Any comprehensive model must ex- 
plain the differences between our blobs and LBGs: 1) a 
comparison of their number density and clustering prop- 
erties suggests that blobs could occu py higher density 
environments (e.g.. lYang et al.ll2010D . 2) the extended 
emission of blobs is typically 10 times brighter, 3) the 
Lya-Ha offset of our blobs is lower. The first and sec- 
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ond differences could arise if the blobs host more mas- 
sive galaxies and higher star formation rates than typi- 
cal of star-forming galaxies at high-z. It remains to be 
tested whether star formation in the galaxies embedded 
in blobs can produce enough Lya photons to explain the 
higher luminosity and surface brightness of Lya halos. 
The Lya-Ha offset is clearly tied in some way to the ex- 
tended Lya emission of our blobs, as it is unusually low 
for both blobs compared to the offse t distribution for 
LBGs (see also iMcLinden e t al. 2010, for compact Lya 
emitters). It would interes t ing to check whether the LBG 
subclasses in lSteidel et al.l (|2011l ) with brighter extended 
emission have smaller Lya-Ha offsets and weaker IS ab- 
sorption at v r~j Okms -1 than the others. For now, we 
note that the RT outflow model discussed in Section[4]ad- 
equately explains the sense and magnitude of both the 
Lya-Ha and the Ha-i nterstellar absorpti on line offsets 
in LAB02, whereas the lSteidel et al.1 ()2010h model would 
require this agreement of the offsets to be a coincidence. 
Of course, we are only able to perform this test using 
one blob at present. As stressed earlier, a larger sam- 
ple of Lya-Ha offsets is required to test whether smaller 
Lya-Ha offsets are generic properties of Lya blobs. 



5. CONCLUSIONS 

Determining the source of Lya emission in extended 
nebulae at high redshift requires that we measure their 
gas kinematics with an optically-thin line such as Ha 
A6563, which, unlike the optically-thick Lya emission, 
is not much altered by radiative transfer effects and is 
more concentrated about the blob's core. The compari- 
son of the Lya and Ha line centers then constrains any 
infall or outflow of the Lya-emitting gas relative to the 
systemic velocity of the blob defined by Ha. We obtain 
optical and NIR spectra of the two brightes t Lya blobs 
CDFS-LAB01 and CDFS-LAB02) from the lYang et all 
2010) sample using the Magellan/MagE optical and 
VLT/SINFONI NIR spectrographs. Both the Lya and 
Ha lines confirm that these blobs lie at the survey red- 
shift, z ~ 2.3 (A c « 4030A). 

The blobs contain several Ha sources, which roughly 
coincide with galaxies in HST rest-frame UV images. 
The Ha detections of these embedded galaxies re- 
veal large internal velocity dispersions (a v = 130 — 
190 kins" 1 ). Furthermore, in the one system (LAB01) 
where we can reliably extract profiles for two different Ha 
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sources, their velocity difference is Av ~ 440 km s -1 . The 
presence of multiple galaxies within the blobs, and those 
galaxies' large velocity disp ersions and large re lative mo- 
tion, supports the result of I Yang et~a l. (2010) that Lya 
blobs inhabit massive dark matter halos that will evolve 
into rich clusters today. The embedded galaxies may 
represent the precursors of brightest cluster galaxies. 

To determine whether the gas near the embedded 
galaxies is predominantly infalling or outflowing, we com- 
pare the Lya and Ha line centers. Lya is not off- 
set (A^Lya = +0kms _1 ) in LAB01 and redshifted by 
only +230kms _1 in LAB02. These offsets are small 
compared to those of Lyman break galaxies, which av- 
erage +450 km s -1 and extend to about +700kms _1 
(jSteidel et alJ 12010). By comparing the observed Lya 
profiles with simple radiative transfer models, we con- 
strain the outflow velocity from the embedded galaxies 
to be ~ - 150 km s _1 in the two blobs. In LAB02, we de- 
tect low-ionization, interstellar absorption lines, includ- 
ing C II A1334 and Si II A1526, whose blueward shifts of 
~ 200 km s -1 are consistent with the small outflow im- 
plied by the redward shift of Lya. Models with simple 
infall geometries or that require outflow velocities ex- 
ce eding ~500 km s -1 (e.g., s uper/hyper-winds discussed 
in iTaniguchi fc Shioval l200"ot) are ruled out. Because of 
the unknown geometry of the gas distribution and the 
possibility of multiple sources of Lya emission embed- 
ded in the blobs, a larger sample and more sophisticated 
models are required to test other infall and outflow sce- 
narios. 
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